Lifetime of the Bose Gas with Resonant Interactions 
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We study the lifetime of a Bose gas at and around unitarity using a Feshbach resonance in 
lithium 7. At unitarity, we measure the temperature dependence of the three-body decay coefficient 
L3. Our data follow a Lz~ ^z/T^ law with A3 = 2.5(3)stat(6)sys x 10~^''(/iK)^cm^s^^ and are in 
good agreement with our analytical result based on the zero-range theory. Varying the scattering 
length a at fixed temperature, we investigate the crossover between the finite-temperature unitary 
region and the previously studied regime where \a\ is smaller than the thermal wavelength. We find 
that L3 is continuous across resonance, and over the whole a < range our data quantitatively agree 
with our calculation. 
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Recent advances in manipulating cold atomic vapors 
have enabled the study of Fermi gases at the unitary 
limit where the scattering length a describing two-body 
interactions becomes infinite. It has been demonstrated 
both experimentally and theoretically that in this limit 
the system is characterized by a scale invariancc leading 
to remarkably simple scaling laws [1]. By contrast, most 
experimental results on Bose-Einstein condensates were 
obtained in the weakly interacting regime. Recent ex- 
perimental results on bosons near Feshbach resonances 
have revived the interest in strongly interacting bosons 
[2]: the development of experimental tools has enabled a 
precise test of the Lee- Huang- Yang corrections [3, 4], and 
several theoretical papers have studied the hypothetical 
unitary Bose gas at zero [5 8] or finite [9] temperature. 

Experimental investigation of ultracold bosons near 
unitarity has been hampered by the fast increase of the 
three-body recombination rate close to a Feshbach reso- 
nance [10, 11]. In this case, the number of trapped atoms 
N{t) follows the usual three-body law 

N = -Ls{n^)N, (1) 

where (n^) = J cPrv?{r)/N is the mean square den- 
sity and L3 is the three-body loss rate constant. In 
the zero-temperature limit Z/3 increases as ha^/m [12] 
multiplied by a dimensionless log-periodic function of 
a revealing Efimov physics [15 24]. At finite tem- 
perature, L3 saturates when a becomes comparable 
to the thermal wavelength Ath = h/ \/2'Kmk-BT, and 
L-i^ha^/m^TV' /m^{kBTf [9, 25, 27]. This satu- 
ration suggests that a non-quantum-degenerate Bose 
gas near a Feshbach resonance will maintain quasi- 
thermal equilibrium [9]. Indeed, in this regime, 
jaj^Ath and nA^j^-Cl. Thus, the elastic collision rate 



^20chXthn/m is much higher than the three-body loss 
rate ^3 = L^v? (xhX^^'n? /m. Experimental and numer- 
ical evidence for a saturation of Lj, were reported in 
[3, 20, 25]. A theoretical upper bound compatible with 
this scaling was derived in [26] assuming that only the 
lowest 3-body hyperspherical harmonic contributes, an 
assumption which breaks down when |a| exceeds Ath- 

In this Letter, we measure the temperature depen- 
dence of the unitary three-body recombination rate and 
find agreement with a L3 a 1/T^ scaling law. In a sec- 
ond set of measurements performed at constant temper- 
ature we study L3 versus a. We show how this function 
smoothly connects to the zero-temperature calculations 
when jaj <C Ath- These observations are explained by a 
general theoretical result for L^{a,T), exact in the zero- 
range approximation, that we derive in the second part. 
Our theory allows for a complete analytic description 
of the unitary case and, in particular, predicts (weak) 
log-periodic oscillations of the quantity L^T^. Our find- 
ings quantify the ratio of good-to-bad collisions in the 
system and provide a solid ground for future studies of 
strongly interacting Bose gases. Furthermore, on the 
a < side, experiments have so far detected a single Efi- 
mov trimcr [3, 21-23, 30]. Our analysis predicts that 
a second Efimov trimer of very large size should be de- 
tectable in ^Li at temperatures on the order of a few 
microkclvins. 

Our experimental setup was presented in [28]. Af- 
ter magneto-optical trapping and evaporation in a loffe 
magnetic trap down to ~ 30 /iK, ~2xl0^ ^Li atoms are 
transferred into a hybrid magnetic and dipole trap in the 
state |1,1)- The transverse confinement is obtained by 
a single laser beam of waist 43(l)/im and wavelength 
1073 nm, while the longitudinal trapping is enhanced 
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FIG. 1: Time dependence of the atom number (a) and tem- 
perature (b) for [/ = r; ksT, with T = 5.2(4) ^K, -q = 7.4 
and (uncorrected) L3 = 1.2(2)stat x 10"^^cm''s"^ The dot- 
ted line shows the long time t~^^^ dependence of the number 
of atoms. 



by a magnetic curvature. The resulting potential has 
a cylindrical symmetry around the propagation axis of 
the laser and is characterized by trapping frequencies 
0.87kHz <a;p/27r< 3.07 kHz and 18 Hz < a;2/27r < 49 Hz. 
Further cooling is achieved by applying a homogeneous 
magnetic field B ~ 718 G for which the scattering length 
is ~ 200ao, and decreasing the depth of the trapping po- 
tential down to a variable value U' allowing us to vary 
the final temperature of the cloud. Afterwards, the dipole 
trap is recompressed to a value U >U' , to prevent sig- 
nificant atom loss due to the enhanced evaporation rate, 
see below. At each T we choose U so as to maintain 
the temperature constant during the three-body loss rate 
measurement. Finally, the magnetic field is ramped in 
100-500 ms to Bo -737.8(3) G where the scattering 
length a diverges [4]. We then measure the total atom 
number N remaining after a variable waiting time t and 
the corresponding T, using in situ imaging of the thermal 
gas. 

Our data are limited to the range of temperature 
< T < lO/zK. For T>l/LtK, the rate -f3 = -N/N re- 
mains small with respect to other characteristic rates 
in our cloud (elastic scattering rate, trapping frequen- 
cies), which guarantees that a quasi-thermal equilib- 
rium is maintained. We check that for these pa- 
rameters the in situ integrated density profile is in- 
deed gaussian, and we use it to extract the temper- 
ature of the cloud, found to be in agreement with 
that of time of flight. The peak phase-space den- 
sity varies within 0.07 x 10"^ < uqX^i^ < 1.1 x 10^^. A 



typical time dependence of N and T is shown in 
Fig. 1. The time dependence of the atom num- 
ber is fitted using the usual three-body recombina- 
tion law Eq. (1) [35]. For a non-degenerate gas of 
temperature T, the density profile is gaussian, and 
we have (n^) ^ N^A{T) ^ N^lmCo'^ /2TrV3kBTf, with 
u) = {ujpUJzy^^ being the mean trapping frequency. We 
then have: 

N = ~L3{T)A{T)N^. (2) 
Assuming constant temperature, integrating Eq. (2) gives 

iV(0) 



N{t) = 



,/l + 2A{T)L,{T)N^{0)t' 



(3) 



which we use as a fitting function to analyze N{t), and 
extract L^iT) as shown in Fig. 1. 

Because of their n^/T"^ dependence, three-body losses 
preferentially remove atoms of low kinetic energy and 
those located at the center of the trap where the density 
is the highest and potential energy is the smallest. As a 
result, three-body loss events heat up the cloud [12]. We 
ensure constant temperature by operating with a typi- 
cal trap depth U ~ rjk-QT with 6 < < 8, for which the 
residual evaporation then balances recombination heat- 
ing, see Fig. lb. This ensures that I/3 is time inde- 
pendent, but, as a drawback, evaporation contributes to 
losses. To quantify the relative importance of evapora- 
tive and three-body losses, we first note that an atom 
expelled by evaporation removes on average an energy 
~ (77-l-K) fceT, where, following [13], we take n from [14], 
while each three-body event leaves on average an excess 
heat of Sk^T per particle. Extending the derivation of 
[12] to the case of an energy dependent three-body loss 
rate cx£'~^, we obtain 5 = 5/3 [36]. The energy bal- 
ance required to keep the temperature constant thus im- 
plies that the evaporation rate is ~5/(ry -t- k — 3) times 
smaller than the three-body loss rate. Neglecting this 
effect would induce a systematic overestimation of of 
about 50% for r] = 6 and 30% for 77 = 8. Therefore, we 
apply this systematic correction to our data. 

The temperature dependence of L3 obtained from 
our measurements at unitarity is shown in Fig. 2. 
It is well fit by the scaling law L^{T) = X3/T^, with 
A3 = 2.5(3)stat X 10~20(^K)2cm6s-i as the best-fit value. 
In order to discuss the systematic uncertainty of this 
measurement we note that the quantity L^T^ scales in 
all experimental parameters identically to the thermody- 
namic quantity {y? / P)"^ of a zero-temperature BEG, with 
chemical potential /i, and pressure P [36]. We use this 
relation to calibrate our experimental parameters [4] and 
obtain a systematic uncertainty on A3 of < 25 % resulting 
in A3 = 2.5(3)stat(6).ys x IQ-^^'iixKYcuY^s-^ 

We now study the a-dependence of ^3 on both sides of 
the resonance by employing the same experimental proce- 
dure as in the unitary case. We tune the scattering length 
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FIG. 2: Temperature dependence of the three-body 
loss rate L3. FiUed circles: experimental data; green 
dashed line: best fit to the data L3{T) = X3/T^ with 
A3 = 2.5(3)stat(6)sys X fO-2''(^iK)2cm«s-i; the green band 
shows the la quadrature sum of uncertainties. Solid line: 
prediction from Eq. (5), A3 = 1.52 x 10^^°(/iK)^cm^s~^ with 
77* =0.21 from [29, 30]. 



while keeping the temperature within 10% of 5.9 /iK, 
see Fig. 3. The excess heat S entering in the correction 
now depends on the value of ka. The correction is ap- 
plied to all data points (filled circles) except in the range 
1500 ao < a < 5000 oo (open circles), where the assump- 
tions of our model are not applicable [36]. In the limit 
|fl|3>Ath, we observe that L^^a) saturates to the same 
value on both sides of the resonance. In the opposite 
limit |a| ^ Ath, our data connect to the zero temperature 
behavior [18] studied experimentally in [20-24]. On the 
a < side, the dashed line is the zero-temperature predic- 
tion for L3 from [18]. We clearly see that finite tempera- 
ture reduces the three-body loss rate. On the a > side, 
temperature effects become negligible for a < 2000 ap as 
testified by our measurements performed on a low tem- 
perature Bose-Einstein condensate (green squares) which 
agree with the total recombination rate to shallow and 
deep dimers calculated at T = in [18] (dashed line). The 
data around unitarity and on the a < side are seen to be 
in excellent agreement with our theory Eq. (4) described 
below. 

In order to understand the dependence L^{a,T) theo- 
retically, we employ the iS-matrix formalism developed 
in [18, 31, 32]. According to the method, at hyper- 
radii i?^ |a| one defines three- atom scattering channels 
(z = 3,4, ...) for which the wavefunction factorizes into a 
normalized hyperangular part, $i(i?), and a linear su- 
perposition of the incoming, R~^/'^e~^''^, and outgoing, 
j^-5/2g+i/ci?,^ hyperradial waves. The channel z = 2 is 
defined for a > and describes the motion of an atom 
relative to a shallow dimer. The recombination or re- 
laxation to deep molecular states (with a size of order 
the van der Waals range Re) requires inclusion of other 
atom-dimer channels. In the zero-range approximation. 
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FIG. 3: a) Scattering-length dependence of the three-body 
rate constant L3(a) for constant r = 5.9(6) nK (filled and 
open circles). For small positive a, L3(a) for a low temper- 
ature condensate is also shown (green squares). The solid 
blue line corresponds to our theoretical prediction Eq. (4) 
for r = 5.9 ^iK. The blue range is the same theory for 5.3 
to 6.5 fiK. The dashed lines represent the zero-temperature 
prediction for L3(a) [18]. All theory curves are for the ^Li 
parameters t;, =0.21 and i?o = 270ao measured in [29, 30]. 
The vertical dotted lines correspond to |a|/Ath = 1. The open 
circles in the range 1500 ao <a< 5000 ao are not corrected for 
residual evaporation as our model is not applicable, b) Loga- 
rithmic plot of the a < side, displaying the two Efimov loss 
resonances. 



valid when i?e ^ = Min (1/k, |a|), the overall effect 
of these channels and all short-range physics in general 
can be taken into account by introducing a single Efimov 
channel (« = 1) defined for Rg <C R<^Rm'. the wavefunc- 
tion at these distances is a linear superposition of the 
incoming, <i>i(E)i?~^+*^'', and outgoing, <i>i(E)i?~^"**", 
Efimov radial waves. Here sq ~ 1.00624. The notion "in- 
coming" or "outgoing" is defined with respect to the 
long-distance region Rm1^R^\<Ai that, for exam- 
ple, the incoming Efimov wave actually propagates to- 
wards larger R whereas incoming waves in all other chan- 
nels propagate towards smaller hyperradii. The ma- 
trix Sij relates the incoming amplitude in the i"^ chan- 
nel with the outgoing one in the j"^ channel and de- 
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scribes the reflection, transmission, and mixing of chan- 
nels in the long-distance region. This matrix is unitary 
and independent of the short-range physics. The short- 
range effects are taken into account by fixing the rela- 
tive phase and amplitude of the incoming and outgoing 
Efimov waves i?^* a (i?/J?o)'*°-e^''* (-R/i?o)"'*°, where 
Rq is the three-body parameter and the short-range in- 
elastic processes are parametrized by r;* > 0, which im- 
plies that the number of triples going towards the re- 
gion of Re is by the factor e^*** larger than the 
number of triples leaving this region [33]. Braaten et 
al. [32] have shown that for a given incoming channel 
i>2 the probability of recombination to deeply bound 
states is P^^{l- e-'^''*)|sii|V|l + {kRo)-^'''°e-^'^' sii\^ 
[34]. For a<0, by using the fact that Sn is unitary 
(Si^i klip — 1) and averaging over the Boltzmann dis- 
tribution we then obtain the total loss rate constant 



72y^TT^h{l-e-''"') f°° (l-|sii|2)e-'='/'=?>.A;(ifc 



'/ 

Jo 



|l-h(A:i?o)-=^"°e-2').sii|2' 



(4) 

where fcth = y/mkBT/h. 

Note that in deriving Eq. (4) we closely followed [32] 
where the scattering length was assumed to be finite. 
However, we easily generalize this derivation to the case 
a = 00, in which the channels become decoupled at dis- 
tances R^l/k and the long-distance region can now be 
defined by i?~ A less trivial result of our analysis is 
that for any ka there exists a unitary transformation of 
the matrix Sij which leaves the clement sn invariant, but 
all channels with i>2> become decoupled from the Efi- 
mov channel [36]. This transformation constructs a new 
large- i? channel characterized by a certain hyperangular 
wavefunction $3(7?). For negative or infinite a this is the 
only channel that can "talk" to the lossy short-distance 
Efimov channel via a unitary 2x2 matrix. Therefore, the 
three-body loss rate can not exceed the so-called maxi- 
mum value Lf^ = SdV^n'^h^ikBT)-'^ /rri^ reached in the 
case when the outgoing flux in this newly constructed 
channel vanishes. Previous derivations of L'^^^ [26] es- 
sentially implied that ^3{R) is the lowest non-interacting 
hypcrsphcrical harmonics. This approximation can be 
made only for fcjaj ^ 1. In general, '^3{R) is not an eigen- 
state of the angular momentum operator. In particular, 
at unitarity l>3(.R) = $i(.R) [36]. 

The function sii(A;a) is calculated in [36]. At unitarity 
it equals suioo) = -e-'^'°e'^'^'°^'''^+^'^^^^+'^°^^ and from 
Eq. (4) one sees that LsT"^ should be a log-periodic fimc- 
tion of T. However, due to the numerically small value 
of )sii]«0.04, in the case of three identical bosons the 
oscillations are very small and Z/3 is well approximated 
by setting Sn = 0: 



36\/3 7r2' 



[kBTf 



(5) 



This explains the L^cxT ^ experimental observation seen 



in Fig. 2 at unitarity. Taking 77* =0.21. which is the 
average of two measurements made for our ^Li Fes- 
hbach resonance in [29, 30], we get L3 = Xs/T'^ with 
A3 = 1.52 X lG-20Q^K)2cm*^s-i. This is 40% below the 
experimentally determined value without any adjustable 
parameter and the agreement between theory and exper- 
iment is 1.4cr. We should point out that Eq. (4) can be 
easily generalized to the case of other three-body systems 
with smaller sq- Then, the terms neglected in Eq. (5) 
can become important. They also become important in 
our system of three identical bosons when departing from 
resonance in the direction of a < 0. Then [sn (fco)| mono- 
tonically increases as a function of l//c]a] reaching 1 in 
the limit /ca— )-0~, the argument of sn also being a mono- 
tonic function of l/fc|a| [36]. The solid blue line in Fig. 3 
is the result obtained from Eq. (4) using the same 77* 
as above and i?o = 270ao also taken from [29, 30]. The 
shaded blue area reflects our experimental range of tem- 
peratures. More or less visible maxima of L3 appear 
when the denominator in the integrand of Eq. (4) reaches 
its minimum, i.e., becomes resonant. The approximate 
condition for this is Argsii(A:a) =7r+2so lnA:i?o a-nd the 
features become increasingly more pronounced for larger 
jsii I and smaller 77*. Note that from the viewpoint of the 
visibility of the maxima, decreasing jaj is equivalent to 
decreasing ^/T. Fig 3b shows the pronounced resonance 
at a = a_ps — 274ao observed in [29, 30]. This resonance is 
associated with the passage of an Efimov trimer through 
the three-atom threshold. Another Efimov trimer, larger 
in size by a factor of e'^/'*° = 22.7, is expected to go 
through the threshold at around a«— 6350ao leading 
to another zero energy resonance. As we deduce from 
Eq. (4) and show in Fig. 3 for 5.9/uK, the thermally av- 
eraged remnants of this predicted resonance lead to a 
maximum of L3 at a~— SlOOao. As seen in Fig. 3b the 
agreement between theory and experiment is very good 
over the entire a < range. 

Because of the existence of a shallow dimer state, the 
case a>0 becomes, in general, a complicated dynami- 
cal problem which should take into account the atom- 
dimer and dimer-dimer relaxation as well as various non- 
universal factors: the finite trap depth, chemical imbal- 
ance between trapped shallow dimcrs and free atoms, and 
deviations from thermal equilibrium, possibly dependent 
on the preparation sequence. These issues require an ex- 
tensive discussion beyond the scope of this Letter. The 
situation obviously simplifies in the case of very small a 
when the system is purely atomic and the three-body re- 
combination to deep and shallow molecules leads to an 
immediate loss of three atoms. 

Discussing the opposite limit of large a>0, we first 
note that dimers are well defined when their size ^ a 
is smaller than n~^^^, which we assume in the follow- 
ing (the limit na^ » 1 is equivalent to the case a = 00). 
In the regime a ^ Ath we find using the Skorniakov-Ter- 
Martirosian equation that S12— ^-O for ka^oo, which im- 
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plies that the atom-dimer relaxation rate vanishes and 
that shallow dimers remain at chemical quasi-equilibrium 
with the decaying atomic ensemble [36]. Shallow dimer 
formation and breakup arc then balanced, so that the 
atomic decay is just given by Eq. (1). The expression of 
Lz for a > was obtained in [32] and reduces to Eq. (4) for 
Si2 0. We conclude that the loss rate must be continu- 
ous across the resonance, in accordance with our exper- 
imental data. Therefore, in Fig. 3a the result of Eq. (4) 
is simply continued to positive a for Ath- 

In summary, we have systematically studied the de- 
pendence of the three-body loss rate on T and a in a 
Bosc gas near unitarity. The T"^ temperature depen- 
dence seems to contradict the well-known Arrhenius law 
Lz oc exp(— AiJ/fesT) for a chemical reaction with ac- 
tivation energy AE. This paradox is resolved by not- 
ing that three-body recombination occurs through the 
— 1/B? Efimov three-body potential that does not display 
any potential hill. Moreover, a quantum interference ef- 
fect in Efimov three-body scattering leads to an enhanced 
decay rate at a negative a, suggesting the possibility to 
observe the signature of a second Efimov trimcr of large 
size. Another future direction is to explore the approach 
to the quantum degenerate regime and test whether the 
virial expansion of the unitary Bose gas can be measured 
by using quasi-equilibrium thermodynamics [9]. 
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I. CALCULATION OF sii(fca) 

The function sii(ka) has been calculated for 
< ka < 10 in [1]. Here we calculate it for ka of any 
sign and magnitude including infinity. The three-body 
Schrodinger equation reads 



[-V^ + K(R)-fc'] *(R) = o, 



(SI) 



where 14 is the sum of binary interaction terms 
parametrized by a and the six-dimensional vector R is 
{(2r3 — ri — r2)/\/3, ri - r2}, where ri,r2 and rg arc 
atomic coordinates, and we set % = m = 1. By defini- 
tion the first row of matrix Sij gives us the solution of 
Eq. (SI) with the following asymptotes. For small R we 
have 

*(R) « $i(^) [(fci?)'*" + s^x{kR)-''°] /s/WqR^. (S2) 

Note that our definition of Sn differs from the one of [1] 

by the factor -{kaf'^'^e''^'^", where (J,, » 1.588 [2]. In 
the asymptotic region of large R we have 

oo 

*(R) « si2*2(R) + siiMR)e''''^/V2kRR^. (S3) 

In Eqs. (S2-S3) all ^i{R) arc symmetrized and normal- 
ized and ^'2(R) is the symmetrized wavefunction of the 
atom-dimer relative outgoing motion normalized to a 
unit flux. Physically ^'(R) describes the stationary flow 
of atoms which arc injected at the origin and can either 
return back with amplitude sn (second term in the right 
hand side of Eq. (S2)) or travel to infinity by using chan- 
nels with i > 3. 

Before explaining the numerical method of calculating 
sii{ka) let us discuss some properties of this function 
which can be derived analytically. In order to do this it is 
convenient to use the complex scaling of the Hamiltonian 
[3] and multiply (rotate in the complex plane) all spatial 
coordinates by the complex number ke~™/'^, i.e., we in- 
troduce R = R/ce^*'^/^. Then the problem reduces to 
calculating properties of the bound trimer state with en- 
ergy E = —1, interaction between the atoms being char- 
acterized by the imaginary scattering length a = ka/i. 
Applying the complex scaling to the asymptotes (S2-S3) 
we see that now the solution is constrained to decay at 
large distances and to have the short-distance asymptote 
oc R'^° + Sii{ka)e'^^° R~^^° . That a is imaginary simpli- 
fies the task: if is a solution for a given real value of 
a, then is the solution for a = —a. This leads to 



the relation sii{—ka) = e~^'^''° /sli{ka), and it is thus 
sufficient to deal, for example, only with ka < 0. Since 

our problem is the inverse to finding the Efimov spec- 
trum versus the three-body parameter and a, the point 
ka = 00 and its vicinity can be treated analytically: the 
wavefunction of an Efimov trimer at unitarity is propor- 



tional to $i(i?) 



which gives 



J,so{^R)-e-^'''J^,soi^R) 

si;^(oo) = -22»«oe-''""r(l+iso)/r(l-iso), the result pre- 
sented in the main text. Moreover, by using the known 
analytic formula for the shift of the trimer energy at small 
1/a with a fixed three-body parameter [4] one obtains 
sii{ka » 1) « sii(oo)(l — Cso/ka), where [4] 



C 



TT sinh 



, Son /SoTT\ 

= 2.1126716... 



47r /so7r\ 

— =cosh 

3^3 V 6 / 



(S4) 



Finally, the three-body wavefunction in the limit of van- 
ishing total energy has been studied in [2] from which we 
obtain sn « {klal)^'"" e''^'^'' in the limit fca 0". 

In order to calculate Si2 for arbitrary ka let us intro- 
duce the reduced wavefunction /(r) defined by 

/ f(2r3 - ri - r2)/y3l = 477 lim |ri-r2|*(R), (S5) 

L J ri->-r2 

and write down the Skorniakov-Tcr-Martirosian (STM) 
equation for the Fourier transform of /(r) (for more de- 
tails see [5]) 



where the integral operator L is defined by 



(S6) 



In 



+p2 - 3fcV4\ 4f{p')p'dp' 

Sk^/AJ 



^/2 



- — PP' 



\/^1Tp 

(S7) 

Note that in Eq. (S6) we use /(p) = f{p) since, in the 
case of three identical bosons, higher spherical harmonics 
of this function correspond to the non-Efimovian kine- 
matics, do not contribute to the asymptote (S2), and do 
not lead to (strong) recombination losses. 

As usual, the branches of the logarithm and of the 
square root are chosen as if the momentum k (or energy 
fc^) is slightly shifted into the upper complex half-plane, 
or, alternatively, p and p' are slightly shifted to the lower 
half-plane. In fact, the complex scaling discussed above 
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means that we rotate p and p' in the clockwise direction 
all the way to the negative imaginary axis and rescale 
them by k. Then changing variables in such a way that 
the integration goes along the positive real axis, we ob- 
tain the same Eq. (S6) in which k ^ i and 1/a — )■ i/ka. 
The resulting equation does not have singularities on 
the real axis and is extremely easy to solve numerically. 
The large-p asymptote of the solution can be written as 
Cip"^"**" + C2P^'^~^^^° , and it is straightforward to show 
that sii(A:a) = (C2/C1) [r(l + iso)/r(l - iso)] e"'^"". 

In Fig. SI we plot |sii| and Argsn versus ka <0 (solid 
line). The dashed and dotted lines correspond, respec- 
tively, to the limits k\a\ 1 and k\a\ <C 1 discussed 
above. 




100 



FIG. SI: Modulus and phase of sn versus ka < 0. Dashed 
and dotted lines are analytic limits, see text. 



Therefore, in order to calculate S12 one has to solve 
Eq. (S6) (with the correct boundary condition at large 
p) and find the residue of this pole. Here we solve this 
problem perturbatively using l/Zca as a small parameter. 

The solution of the three-body problem at unitarity is 
given in terms of the BcsscI functions. The corresponding 
correctly normalized function fo{p) up to a phase factor 
can be written as 



27r-\/Cso cosh(7rso) exp(— ttsq) 



(SIO) 



\/3 sinh(7rso /2)p^/lP - k'^ 



In deriving Eq. (SIO) we used the small R asymptote of 
the normalized hyperangular Efimov wavefunction $i(i?) 
which we took from [4]. The constant C is defineed in 
Eq. (S4). 

Let us now write the solution of Eq. (S6) at small 1/ka 

f{p) = foip) + ^f{p): where 5f{p) tends to zero when 
a — > 00. Equation (S6) now reduces to 



fc2 - l/a)Sfip) = Up)/ a + Lk^dfip). (SIX) 



II. LARGE ka ASYMPTOTE OF \si2{ka)\ 



Looking at the right hand side of this equation at p = po 
we observe that the first term tends to a finite value as 

a — > 00 since /o(p) is singular at p — > fc. In contrast, 
the integral operator smooths singularities and makes the 
second term vanish uniformly for large a. Therefore, the 
dominant contribution to S12 can be obtained by neglect- 
ing the second term, and we finally obtain 

|si2(fca)|2 « 2Cso [1 + exp(-27T.so)] /ka, ka > 1, (S12) 

i.e., we have managed to show that S12 — >■ 0, as one 
approaches the resonance. 



The quantity ,si2 is the amplitude of the atom-dimer 
outgoing wave, see Eq. (S3) . By construction it is zero for 
a < since there are no shallow dimers and, as we argue 
in the main text and in Sec. Ill, for a > it becomes 
important for determination of L3. The question that 
we address now is whether |si2p vanishes for a — )■ -|-oo 
or not. Unfortunately, numerical results for Si2{ka) are 
available only for ka ^ 10 [1] and do not allow us to make 
any statement on the large ka behavior of this quantity. 

In the region |ri — r2| ~ a and R ^ a the atom- 
dimer wavefunction ^'2(R) introduced in Eq. (S3) can 
be written as 



*2(R) 



exp(-|r2 - ri|/a) exp{ipoR) 
\/27ra|r2 - rJ ^/24^^p^R' 



(S8) 



where po = \/k'^ + > fc is the atom-dimer relative 
momentum. The outgoing wave (S8) corresponds to the 
pole of f(p) at p = Pi:,: 



2TTSi2{ka) 



V^Sopf P-Po 



\p-Po\^Po-k. (S9) 
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III. ATOM-DIMER CHEMICAL EQUILIBRIUM 
NEAR RESONANCE 

On the positive side of the resonance the loss rate in 
the system is no longer solely due to the recombination to 
deep molecular states. Three atoms can recombine to a 
shallow dimer and depending on how its binding energy, 
Ed = h^/ma'^, compares to the trap depth, U, the prod- 
ucts of such a three-body event may or may not leave 
the trap. Moreover, even if they have enough energy to 
leave, they can collide with the remaining atoms and re- 
distribute their excess energy into heat. This dynamical 
problem, in general, goes far beyond calculating the loss 
rate in three-atom collisions. Obviously, this complica- 
tion is absent very far from the resonance where, starting 
from a purely atomic sample, one counts any recombina- 
tion event as the loss of three atoms. By contrast, if 
Ed <U , the shallow dimers stay in the system and mix 
with atoms. 

We now focus on the regime Ed <C fc^T, where the 
situation greatly simplifies. Let us assume chemical equi- 
librium between atoms and shallow dimers, and validate 
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this assumption a posteriori. The dimer density n£) is 
then related to the atomic density n by 

no = n'Xl^ 2V2 e^^/^'^-^) n^Xl^ 2^2. (S13) 

The two reverse processes of three-atom recombination 
to a shallow dimer and of atom-dimer breakup (i.e. dis- 
sociation of a shallow dimer after collision with an atom) 
then balance each other and give a vanishing total con- 
tribution to dn/dt. Hence dn/dt = —L3 — nriD, 
where L3 is the rate constant for recombination to deep 



dimers and L2 is the rate constant for atom-dimer re- 
laxation (i.e. formation of a deeply bound dimer after 
collision of a shallow dimer with an atom). The expres- 
sion of Ls{T) for a > was obtained in [1] and differs 
from Eq. (4) of the main text only by the replacement of 
the term 1 — |siip by 1 — |siip — |si2p- For the incom- 
ing atom-dimer channel i = 2, the expression of the loss 
probability P2 was given above Eq. (4) in the main text, 
and leads after thermal averaging to 



^2 ^^>~ {mkBTf/^ 



S12 



|l + (|fc|i?o)-2-oe-2'7.sn| 



(S14) 



Here the integration variable E = fi^k^/m is the total 
energy of the three-atom system in the center of mass ref- 
erence frame. Thus, the integration over negative E de- 
cribes the atom-dimer relaxation events below the break- 
up threshold. We should also note that the matrix ele- 
ments sii and S12 are functions of VSa, which becomes 
imaginary for < 0. To show that vanishes in the 
large a limit, we treat separately the contributions from 
positive and negative E: For > 0, we have seen above 
that S12 vanishes as 1/a, see Eq. (S12), which leads to 
a contribution cx 1/a to i^"^; for E < 0, the integra- 
tion is limited to the narrow window [—E]j,0\, and the 
integrand can be bounded from above thanks to P2 < 1, 
leading to a contribution oc 1/a^ to ^2^. This allows 
us to neglect in the rate equation for dn/dt, which 
then reduces to Eq. (1) of the main text [14]. 

Finally, let us validate our chemical equilibrium 
assumption. For a given dimer, the event rates 
for relaxation and breakup (after collision with an 
atom) are respectively 7rci = nL2^{T) and 7broak = 
nashaiiow(r) y2e-^°/('=^^), where ashaiiow(T) is the 
event rate constant for three-atom recombination to a 
shallow dimer, and we used Eq. (S13). In the regime 
a ^ Ath considered here, we can estimate from [1] that 
3ashaiiow(T) saturates to a value > lOLf^iT) [15]. 
Evaluating the leading-order behavior of L^'^ as ex- 
plained above then gives 7rei/7brcak 0.1 Ath/a ^ 1 [16]. 
Hence the relaxation events do not destroy chemical equi- 
librium, as they happen much less frequently than the 
breakup events (and thus also than the reverse dinier- 
formation events). The rclax;ation rates in dimcr-dimer 
and dimer-atom-atom collisions are also smaller than 
7break by factors oc nAjjj <C 1 (with unknown prefac- 
tors which depend on the four-body problem). The last 
condition to check is 7break ^ 73, i-e., a given dimer 
should be likely to break up (and to be replaced by a 
newly formed dimer) within a time miich smaller than 
the timescale I/73 over which the cloud decays. Esti- 



mating 7break as above and using our result for 73 gives 

Is/lbreak ^ 0.1 TlAf^ < 1- 



IV. MATRIX Sij-. NEGATIVE a 

In the case a < the atom-dimer channel is closed 
and the structure of matrix Sij is as follows. We have a 
single discrete small- i? Efimov channel and a continuum 
of large- i? channels, {$i(i?), i > 3} being a complete 
orthonormal set of hyperangular functions, for example, 
eigenfunctions of the hyperangular kinetic energy opera- 
tor in the absence of interactions. Given Sij we change 
this basis in favor of another orthonormal set {^i{R)} in 
which we choose 

i>3 

The normalization condition which we impose on $3 
uniquely defines this function and S13 (up to an ir- 
relevant phase factor). The corresponding asymptotic 
triatomic channel is defined by the incoming, -03 = 
^3{R)e-'''^/V2kRR'^ and outgoing, VI , waves. Note 
that we do not touch the Efimov channel and, therefore, 
sii = sii. From Eq. (S15) and the imitarity of .s^ we 
can deduce that |sii|^ -|- [sisj^ — 1. Since the new ma- 
trix Sij should also be unitary, we conclude that sn and 
Sis are the only non-zero entries of its first row. Let us 
now write explicitly the corresponding wavefunction [cf. 
Eqs. (S2-S3)] 

where we denote the incoming Efimov wave as Vi = 

^i{R){kRy">/^^R'^. 

Choosing an appropriate linear combination of ^E* and 



4 



loss 




short-distance 
region 



scale-invariant 
region 



long-distance 
region 



asymptotic 
region 



FIG. S2: A thrcc-body wave arriving from large hyperradius 
R with amplitude ^3° in the triatomic channel i = 3 can follow 
various pathways before it either returns to large R, or gets 
lost at i? ~ _Re by turning into an atom and a deep dimer. 
One can imagine a Fabry-Pcrot interferometer, the mirrors 
of which are formed by the short- distance and long-distance 
regions. Multiple reflections by these regions can lead to the 
resonant denominator in the three-body loss rate formula. 



^* we obtain the relation 



'13 



S131PI, R 

- (Sl3Sn/st3)V'3> -R- 



00, 



(S17) 



the right hand side of which defines the second {i = 3) 
row of the matrix Sij. Namely, S31 = S13, S33 = 
-{siasti/sls), and hj = for j > 3. 

Because s is unitary one sees that sn = Sis = for i > 
3. Therefore, the upper left 2x2 block completely sepa- 
rates from the rest of the matrix. The problem of calcu- 
lating the three-body loss rate then reduces to the prob- 
lem of finding the four amplitudes and ^3"''°"', 
which are the coefficients in front of the corresponding in- 
coming and outgoing waves in the three-body wavefunc- 
tion, see Fig. S2. The coefficient ^3" is found by project- 
ing the initial correctly normalized six-dimensional plane 
wave into the state ips. The amplitudes are re- 
lated by the three-body contact condition A^-^ = AA°^^ 
where A ~ — (fci?o)^^*^"e^^''* . Finally, the relation be- 
tween the incoming and outgoing amplitudes given by 
the matrix Sij provides the last two linear equations nec- 
essary to solve the problem: = s\\Af + and 
= S3i^'" + sasAjj". The loss rate is then obtained 
by calculating the difference between incoming and out- 
going fluxes either for i? — > oc' or i? — > 0. Averaging over 
the thermal distribution one recovers the formula for L3 
presented in Eq. (4) of the main text, where we now see 
that the non-trivial fc-dependence of the integrand comes 
from the interference between the various pathways rep- 
resented in Fig. S2. The simplified approximate formula 
given in Eq. (5) of the main text corresponds to neglect- 
ing any reflection from the long-distance region, hence 
no more interferences and no log-periodic modulation of 

^3^2 with Ath/it!t. 



MATRIX Sii AT UNITARITY 



At unitarity, the procedure of transforming the matrix 
Sij into block diagonal form is very simple. Having an 
inflnite a does not introduce a lengthscale into the prob- 
lem and, as a consequence, the adiabatic hyperangular 
eigenf unctions do not depend on the hyperradius, lead- 
ing to the complete separability [7] between the hyperan- 
gular and hyperradial problems. Namely, the three-body 
wavefunction can be written as 



*(R) = ^0,(i?)^^,(i?)i?- 



(S18) 



where (j>s{R) and are, respectively, the (normalized) 
eigenfunctions and eigenvalues of the hyperangular ki- 
netic energy operator supplemented with the unitary 
two-body contact conditions. The hyperradial wavefunc- 
tions satisfy 



dB? 



1 d 
RdR 



s 



F,{R) = eF,{R). (S19) 



In the case of three identical bosons considered here, the 
set {s} contains a single imaginary number s = i Sq — 
il. 00624 (Efimovian sector) and an infinite number of 
real numbers (non- Efimovian sectors). In the Efimovian 
sector the attractive —Sq/R^ potential gives rise to the 
following asymptotic behavior of J^j^o- for R ^ ^/k 
we have Fisg{R) oc R"^^^" and in the opposite limit 
Fisg{R) oc exp{±ikR) / \/R. These two asymptotes of the 
same function actually define the two channels i = 1 and 
i = 3 discussed in Sec. I, the hyperangular wavefunc- 
tions being $1 = $3 = ^jg^j. The rest of (j)s, appropri- 
ately relabelled, form the rest of the set The corre- 
sponding matrix Sij has a 2 x 2 block in its upper left 
corner, which describes the transmission and reflection 
of the wavefunction Fisg (R) by the long-distance region 
R ~ 1/k. The rest of is simply diagonal because (i) 
these channels are decoupled from each other and (ii) the 
repulsive s^/i?^ potentials do not allow (in the zero-range 
approximation) for a transmission of the corresponding 
waves to the short- distance region R ^ R^. Solutions of 
Eq. (S19) can be written in terms of Bessel functions. In 
particular, the wave that has properties of Eq. (S16) can 
be written by setting 



PiSQ (R) 



2i«or(l + iso) 



[J^soikR) - e 



,{kR)] 



(S20) 

in Eq. (S18). Expanding Eq. (S20) at smaU R 
we obtain the result for sii(oo), which has already 
been mentioned. For completeness, from the la rge- 

asymptotes of J±is„ we get .§13 = 2'*° r(l -|- 

iso) sinh(7rso) exp(— 7rso/2 — in/A). 
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VI. MATRIX Sij-. POSITIVE a 

In the case o > we have to take into account another 
discrete channel: the large- i? atom-dimcr one denoted by 
i = 2. By using a similar construction as in the case of 
negative a, one can show that the matrix Sjj can be re- 
duced to a block-diagonal form with a 4 x 4 block in the 
upper left corner, i.e., there are actually two triatomic 
channels, i = 3 and i = 4, coupled to the atom-dimer 
and Efimov ones, and decoupled from the rest of the tri- 
atomic continuum, i > 4. As we have shown in Sec. II, 
Si2{ka) — >■ as a ^ oo. It is then straightforward to 
show that in this limit channels 1 and 3 approximately 
decouple from channels 2 and 4. This means that dimers 



existing in the system are more likely to break-up or scat- 
ter elastically [17] than to relax to deeply bound states. 
This is consistent with our earlier conclusion on the atom- 
dimer chemical quasi-equilibrium close to the resonance. 



VII. EVAPORATION AND 
ANTI-EVAPORATION 

For a<0, the three-body recombination to 
deeply bound states gives the contribution 
Nobody = — j Lzn^{r)d^r to the atomic decay, and 
the corresponding energy loss rate equals 



Sbody 



I 



df^r 



3C/(r) + 



+ 



(r) 72V3TT^h{l - 6-4"*) f°° n'^k'^ (1 - |sii|2)e-'='/'=' fcdfc 



mk^ 



I 



m |l-|-(fci?o)-2'^°e-2'7.sii|2 j 
(S211 



where Lan"^ (r)/3 is the frequency of three-body events 
per unit volume, 3C/(r) and 3A:bT/2 are the loss of 
trapping potential energy and of center-of-mass ki- 
netic energy by each recombining triple, and the 
last term is the loss of relative-motion kinetic en- 
ergy [13]. Let us write the lost energy per lost atom 
as Ezbody / Nobody = (3 - 5)kBT where dksT is the ex- 
cess energy as compared to the average energy per atom 
U:bT. 

For evaporation, Eevap/ Ngvap ~ {1] + i^)kBT, where we 
can take the expression of k in terms of rj given for a 
harmonic trap in [8]. Indeed, as realized in [9], two-body 
collisions leading to an evaporative loss occur mainly in 
the cloud center where the trap is harmonic, and the 
relative momentum for such a collision is approximately 
fixed by the trap depth so that the result derived in [8] for 
an energy- independent two-body cross-section is applica- 
ble. The condition of constant temperature means that 

Eevap + Embody = 3kBT{Nevap + Nobody), which yiclds 

Nevapl Nobody — S/{r] + K — 3). At unitarity, we can ne- 
glect fill in (S21), which gives 6 « 5/3. For — a<ICAth, 
we recover S = l as in [11]. 

For a > 0, we use S = 5/3 when a > Ath- In the opposite 
limit a<C Ath we use 6=1 from [11]. 



ibrate our measurement of the value of A3 [12]. We pro- 
duce a condensate at a = 200 uq and measure the normal- 
ized pressure h versus the gas parameter u = ^a^, where 

g = Airfi^a/m. Next, we find that in order to match h{v) 
to the mean-field prediction, we must multiply the pres- 
sure by a constant ^ = 2.4. ^ corrects for errors in the 
calibrations of our experimental system, e.g. our absolute 
atom coimting, through the product Ur"^ /{ujz'^ {pxf' ao), 
where px is the size of a camera pixel magnified through 
the imaging system to the gas location and (Tq is the 
atomic absorption cross-section for the imaging light. 
When fitting Eq. (3) to our data, we extract two fit 
parameters: 73 = A{T)Lz{T)N'^{id) and 7V(0). Conse- 
quently, our result for A3 = L^T'^ oc -f^T^ /{N^{0)ojr^ojz'^) 
scales as ujz^ (pxf gq^ /ojp'^. This factor is exactly ^~2. 
Ultimately, we estimate our uncertainty in to be 
25 %, dominated by the uncertainty in our trap frequency 
measurements. For the data in Fig. 3, we have an addi- 
tional 20 % uncertainty arising from the scaling with 



VIII. DISCUSSION OF UNCERTAINTIES 

We make use of the grand-canonical equation of state 
for a degenerate Bose gas in the mean-field limit to cal- 
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